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Abstract: We describe string frame bouncing universe scenarios involving the creation 
and annihilation of a non-BPS D9-brane in type IIA superstring theory. We find several 
classes of solutions, in which the bounce is driven by the tachyon dynamics of the non-BPS 
brane. The metric and the dilaton are consistently described in terms of the lowest order 
effective action. The bounce solutions interpolate between contracting and expanding pre- 
big bang (or post-big bang) solutions. The singular behavior of our bounce solutions is 
the same as that of the pre-big bang (or post-big bang) solution. Upon adding a simple 
dilaton potential the asymptotic curvature singularity is removed but the dilaton still grows 
without bound. Such a potential may result from a' corrections in the open string sector. 
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1. Introduction 

The resolution of the big-bang singularity is an important open problem in standard cos- 
mology (see e. g. ||l], ^ for a review and references therein). At the same time it is a natural 
playground for string theory since quantum gravity corrections are expected to be relevant 
in this regime. In the ekpyrotic scenario ^ and a refined version, the cyclic universe 

the hot big bang is the result of the collision of two branes. Explicit cyclic models 
have been suggested as effective four- dimensional models inspired from heterotic M-theory. 
The bounce in the ekpyrotic scenario occurs at a real curvature singularity and thus does 
not resolve the curvature singularity. The new ekpyrotic scenario realizes an explicit 
bounce dynamics by addition of a ghost condensate but suffers from a vacuum instability 
problem (see 0). More generally, due to the necessary violation of the null energy condition 
(NEC) during a bounce, phenomenological models producing a bounce often suffer from 
the problem of introducing matter with negative energy density, i. e. ghosts. 

There are other ways to address the big-bang singularity problem in string theory^ . The 
pre-big bang scenario (see ||ll| for a review) is a consequence of the fact that the tree-level 
equations of motion of string theory are not only symmetric under time reflection t i— > — t 
but also symmetric under the scale-factor duality transformation cl ^ with an appro- 
priate transformation of the dilaton. The 'post-big bang' solution of standard cosmology 
with decelerated expansion defined for positive times is by these dualities connected to an 
inflationary 'pre-big bang' solution for negative times. In this way the cosmic evolution is 



^Alternative suggestions avoiding the problem of introducing matter with negative energy density include 
loop quantum cosmology |^ and matrix models (e. g. |^ [ic||). 



extended to times prior to the big bang in a self-dual way but the solution is still singular. 
One can obtain regular self-dual solutions by tuning a suitable potential for the dilaton. 
But albeit a potential of this form might be the result of higher-loop quantum corrections, 
its form has not been derived from string theory. 

In this paper we consider a novel scenario in string theory where a bounce occurs 
in the string frame due to the creation of an unstable (non-BPS) brane as the universe 
bounces through a string size regime before expanding as the brane decays. Our solutions 
effectively interpolate either between a contracting and an expanding pre-big bang solution 
or between a contracting and an expanding post-big bang solution. The future (past) 
singularity of the pre (post)-big bang solution is not resolved in this scenario. The nice 
feature in our model is that the curvature as well as the dilaton and its derivative remain 
small (in string units) through the bounce so that referring to perturbative string theory 
and the simplest low energy effective action for these fields is justified. In addition no fine 
tuning is required. On the other hand, we do not address issues like dilaton and moduli 
stabilization which are of course important problems to embed this model into late time 
standard cosmology but are not relevant during the string scale regime where the bounce 
occurs. We should stress that our string frame bounce solutions describe monotonously 
contracting or expanding geometries in the Einstein frame. A bounce in the Einstein frame 
may occur upon stabilizing the dilaton asymptotically. However, this entails violating the 
NEC or, alternatively, allowing the gravitational coupling to change sign in the string frame. 
We will discuss a model where the latter effect occurs. 

Our work is organized as follows. This paper consists of seven sections, of which this 
introduction is the first. The model we employ is described in Section ||. In this section 
we present the effective actions and the equations of motions for the metric, dilaton and 
tachyon from a non-BPS brane, which makes the bounce possible. In Section |3| we show 
that our model has no ghost by considering the null energy condition in the Einstein frame, 
and some features of string frame bounce scenarios are studied in relation to the Einstein 
frame. In Section ^, the asymptotic behavior of the solutions is analyzed and its qualitative 
similarity to pre-big bang scenario is clarified. The numerical determination of the global 
bounce solution is presented in Section ^. In Section ^ we present a simple model that 
resolves the asymptotic curvature singularity in the string frame. Finally we conclude in 
the last section. 



2. The model 

We consider a non-BPS space-filling D9-brane in type IIA superstring theory. The details 
of the compactification will not play a role here. Concretely we consider the lowest order 
effective action for the metric and dilaton in the string frame as well as an effective action 
for the open tachyonic mode of the non-BPS D-brane. For now, the only assumption being 
made for the tachyon action is that only the first derivatives of the tachyon appear in it. The 
ansatz for the gravitational action is justified provided the dilaton and metric are slowly 
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varying in string units. We write 

S = -\- [ di°xv^e-2*(i? + 45^$5^$) + 5T (2.1) 

with 5t = j d^°x^e-'^L{T,d^Tdf'T), (2.2) 

where $ is the dilaton, T is the tachyon, and k^q = SvrGio with Gio the ten- dimensional 
Newton constant. We use the signature (— , +,...,+) for the metric. With these conven- 
tions, the matter energy-momentum tensor is given by 

For a Lagrangian minimally coupled to gravity, the metric appears only in (9T)^ = dpTd^T, 
and we can write the energy-momentum tensor as 

= 2 e"* ^^qIIqtY)^ d^Td^T - 5^ e'^L (2.4) 

In a homogeneous isotropic universe that we will consider, all fields are assumed to depend 
only on time; the energy density e and pressure p are given by 

e = To" , p = -Ti (no sum) = e'^L. (2.5) 

We are now ready to make an ansatz for the metric. We take a four- dimensional 
spatially flat FRW spacetime times a six-torus characterized by a single modulus a. Namely 

ds^ = -dt^ + a{tf 5ijdx'dx^ + e^'^^ 5udx^dx\ (2.6) 

where the lower-case Latin indices run over the three uncompactified space coordinates, 
while the upper-case Latin indices label the six compactified dimensions. We further sim- 
plify the problem by restricting to the case where a{t) = e°^^*^. Although equality of the 
two scale factors is phenomenologically not satisfying at late times it may be assumed near 
the cosmological bounce which is the prime focus of this paper. In particular, we will not 
address the important problem of moduli- and dilaton stabilization required to connect to 
the standard cosmology at late times. With the ansatz ( ^ ) the Einstein equations are now 
effectively isotropic. And in particular the relations (^) can be used. From now on the 
indices i, j include I, J. 

The equations of motion for gg, gl and $ are then given by the following first three 
equations 

72H'^ - 36H^ + _ 2k?o e^* e = 0, (2.7) 
2^-8H + 16H^ - 2^2 - 36H^ - kIq e^* p = 0, (2.8) 

2<^> + 18H^ - 2$2 -9H - 45i?2 - ^ e^* p = 0, (2.9) 

€ + 9H{e + p) --l^p = 0, (2.10) 
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and the last equation follows from the generalized conservation law V ^Ty = (9jy<I>) e * L. 
In the case where e = and p = 0, these equations allow exact solutions, namely 



H 



-1±3 



These are the special cases of the pre -big bang {t < 0) and post-big 

-t) and scale-factor 



'-St' 2t 

bang [t > 0) solutions, which respect the time reflection symmetry {t 



duality symmetry (a 



(see for example [1 



Let us explore the possibility of the bounce. Subtracting Eq. ( |2.9D from Eq. ( |2.8| ) we 



find 



H + 9H'^ -2H^- ^e^'^p = 0. (2.11) 
This is an important equation. It tells us that a necessary condition to have a bounce, 



H = 



and 



H>0. 



(2.12) 



is that the tachyon pressure p must be positive. This is a tight constraint for a scalar 
field action; the Born-Infeld action for instance, which is an often used ansatz as a higher 
derivative scalar field action, gives a pressure that is always negative. 

Furthermore, assuming that H is negative during a contracting phase with growing 
dilaton and a negative pressure, ( 2.11j ) implies H < 0, i.e. accelerated contraction. On 
the other hand, for positive equation of state for the scalar field, w > 0, the conservation 



equation in ( 2.10 ) gives a growing pressure p in the contracting phase so that a "turn around" 
iif = is compatible with ( 2.11j ). 

Of course, a Born-Infeld action for the tachyonic sector of non-BPS branes is not in any 
way suggested by string theory. On the other hand, within the restriction to first derivative 
actions it is possible to derive an approximate effective action from string theory for the 



open string tachyon of an unstable brane. This action, constructed in |13| and further 



studied in |14| is given by 



L = -V2Tge-2^ (e-^^^ + \/7r(5r)2 erf ( yW) 



(2.13) 



where rg is the tension of a BPS 9-brane, and therefore V^Tg is the tension of a non-BPS 



9-brane ||15|. Let us shortly summarize how this action was constructed. First, setting 
(5T)2 to zero, we see that the potential is given by 



V{T) = V2Tge 



(2.14) 



This is the exact potential for the open string tachyon potential found in boundary su- 
perstring field theory 18]. The locations of the minima of V{T) are at T = ±00. 

At these values the energy is degenerate with the closed string vacuum which means that 
the non-BPS brane is absent. The construction of the full action ( 2.13| ) is based on the 
observation that the tachyon kink T{x) = x sin(x/\/2a'), where x is one of the spatial 
world volume coordinates, is an exactly marginal deformation of the underlying boundary 



conformal field theory ||19|, |20| and thus should be a solution of the equations of motion 
obtained from ( |2.13 ). It turns out that this requirement determines uniquely the action 
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once the potential has been chosen. Furthermore, it follows by analytic continuation that 
Sen's rolling tachyon solution ||2l|| 

T{t) = A sinh{t/V2a') + B cosh(t/\/2c7) 

is also a solution of the action ( |2.13| ) for all values of A and B. In this dynamical decay 
(or creation) of the non-BPS brane the energy is conserved. The asymptotic state for large 
positive (or negative) times has been argued to be given by "tachyon matter" - essentially 
cold dust made from very massive closed string states. 




Figure 1: X-dependence of the pressure Figure 2: Equation of state w of the 
F{X). tachyon as a function of X. 

Let us now briefly explain how this action can allow a positive pressure, or equivalently 
a positive Lagrangian This follows from the fact that it is real and continuous also 

for negative values of {dT)"^ . Indeed if we write — (5T)^ = X (note that X = in the 
homogeneous case), we can see that 

V^ttX erf (^^) = -2 \/X / e"' ds . (2.15) 

Jo 

This is negative and grows in absolute value faster than the first term e'''^ in the Lagrangian; 
so for positive enough X (for negative enough (9T)^), the Lagrangian, and thus the pressure, 
is always positive. This can be seen from Fig. |l|, where we show the X - dependence of 
the pressure, F{X) = — ^e-^ — 2\fX J(f^ e^^ ds^ . In terms of F the Lagrangian can be 

expressed as L = \/^Tge^2^F{X). From the figure it is clear that dF/dX > (this is 
also clear from ^ = ^ e^^ ds > 0). From now on we work in the unit system with 

a' = 1/2. The energy density and pressure of the tachyon are 

e = V2r9e-*e-^'+^' , p = V2Tge~''e'^" F{f^) . (2.16) 

Note that the equation of state w = e/p depends only on as shown in Fig. ^. In 
particular, — > as ^ oo, while w — > — 1 as — > 0. 

3. Einstein frame and null energy condition 

In this section we show that our model has no ghost and satisfies NEC in the Einstein 
frame. This is important because there might be a pathology due to an instability coming 



- 5 - 



from NEC violation. On the other hand, the consideration on NEC will help us draw more 
general conclusions concerning the bounce scenarios. 

The action (^j|)- (|2.2|) in the string frame with L given in ( ^.13|) is expressed in the 
Einstein frame by means of a conformal transformation g^y = g^,^ e 2 (where 5^,^ is the 
metric in the Einstein frame). This yields 
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2 



d^°xy^e-L(T,Xe-2), (3.2) 

where R and V are the scalar curvature and the covariant derivative associated with g, and 
X = -{VTf = Xet. 

The relevant quantity for verifying classical stability and nonexistence of ghost is the 
sign of the slope of the kinetic term with respect to the first derivative of the fields, which 
should be positive. If it is negative, it signals the existence of a ghost (quantum mechanical 
vacuum instability) and also that the squared speed of sound is negative (classical instabil- 
ity) 0. Since the dilaton has the correct sign for the kinetic term in the Einstein frame, 
the only possible source of ghost is the tachyon Lagrangian. Therefore, one should only 

check the sign of 4^. We have that 4^ = ^44~et^ = ^ /n'^e^' ds > 0. So there 

dX dX Q-^ dX vX '"J 

is no ghost and therefore no violation of NEC in our model. 

Now we turn to the issue of the bounce. In the spatially flat FRW spacetime we are 
considering, the Hubble parameter He in the Einstein frame is related to that in the string 
frame by He = H — ^/A. Since the null energy condition is not violated in our model, 
the Hubble parameter in the Einstein frame monotonously decreases, so that the bounce 
cannot arise in the Einstein frame. If the dilaton were frozen, it would be impossible to 
have a bounce in the string frame as well because the two frames would be trivially related 
(in particular He = H). This is why the running dilaton is crucial for the string frame 
bounce. Now we make some remarks on the bounce scenario in the string frame. We 
assume that the bounce arises as an interpolation between two out of four different phases 
(i. e. contracting or expanding pre/post-big bang phases) in the pre-big bang scenario. Four 
transitions are then possible, namely from the contracting pre-big bang phase to expanding 
pre/post-big bang, or from contracting post-big bang to expanding post /pre-big bang. But 
the transition from pre-big bang to post-big bang cannot happen in our model because the 
NEC in the corresponding Einstein frame is not violated. To show this, we note that in the 
pre-big bang scenario with 5t = the solutions are 

n ■ 9n — 1 1 , , 

H= , ^ = —, r with n = ±-. 3.3 

t-to 2(t-to) 3 ^ ' 

Here, t < Iq corresponds to the pre-big bang phase, and t > to to the post-big bang 
phase. It then follows that He = g^^~^^^ is negative for pre-big bang solutions and positive 
for post-big bang solutions, meaning that pre/post-big bang solutions correspond to a 
contracting/expanding universe in the Einstein frame, regardless of whether the universe 
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is contracting or expanding in the string frame. This means that a transition from pre-big 
bang to post-big bang in the string frame corresponds to a bounce in the Einstein frame, 
which is impossible unless the NEC is violated. Therefore, under the assumption mentioned 
above, the only possible string frame bounce scenarios in our model are the interpolations 
either between two pre-big bang phases or between two post-big bang phases. 

This argument can be extended to the case where the asymptotic behavior of the 
solutions in the string frame is qualitatively, but not exactly, in agreement with the pre- 
big bang scenario. This is the case in our model, as we will see in the next section. In 
conclusion, if one is given a model that does not violate the NEC and if one knows the 
asymptotic boundary conditions of the solution in the string frame, one can then predict 
the possible bouncing scenarios in the string frame by looking at the corresponding Einstein 
frame. An example will be given in the next section. 



4. Asymptotic analysis 

In this section we will try to obtain approximate analytic solutions. This will provide us 
with the asymptotic boundary conditions for the numerical solutions in the next section. 
We emphasize here that by "asymptotic" we mean t — > — oo, or t approaching a pole to, at 
which H diverges, as in the pre-big bang scenario. Similarly the present analysis applies to 
i — > oo, or i approaching a pole to, at which H diverges as in a post-big bang scenario 

To simplify our analysis, we will assume that in either of these limits, the tachyon 
behaves like dust, i.e. p = w{t)e with w{t) 0. This is equivalent to claiming that 



r| — > oo asymptotically because p oc e/(T) when |r| — > oo (see |14| for details). To justify 



this assumption, we look at the tachyon equation of motion following from (|2.13|) . We have 

f + (dH D{f) -T = 0, (4.1) 



where the function D{y) = e~^^ e'^^ds is known as the Dawson integral. This function 
is an odd function, and thus vanishes at y = 0. We will also use the fact that D{y) = 
^ + 0{y~^) for |y| —>■ oo. For |t| — > oo, we will see that {9H — ^) tends to zero. We 
can thus ignore the second term in the equation of motion^ (13) ■ Thus T = T and T, 
as well as T, will grow exponentially; and the pressure will thus vanish for |t| — > oo. We 
emphasize that since the pressure vanishes exponentially fast, we can simply remove it from 
the asymptotic equations of motion because it will always be dominated by the other terms 
that will vanish with a power law. 

When t approaches a pole to the analysis is slightly different because there the term 
{9H — <I>) diverges (unless ^ ~ 9H, but we will see that this does not happen on our 
numerical solution) . We will further assume that this term diverges at least as fast as . 
Assuming that T is finite at to this then implies that either T or T will diverge like or 
faster, in contradiction with a finite T. We cannot exclude the case where r(to) vanishes 
in such a way that it precisely cancels the divergence in {9H — $). In that case the second 



^Note that we ignore the possibility T ^ in the infinite past or future; we are only interested in the 
cases where the D-brane is absent in these limits. 
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term of the equation of motion could be regular at to, and thus T and T could be regular 
there as well. We will nevertheless ignore this possibility because T{to) corresponds to 
very particular initial conditions. For generic initial conditions we will therefore have that 
|r| oo when t —>■ to- This then justifies our claim that we can ignore the pressure for the 



asymptotic analysis. An immediate consequence of Eq. ( 2.10 ) is then that e ~ a~ . 

We now will proceed by analyzing the system of equations (2.7- 2.101 ) assuming that 
either ^ (x H, |$| <C |^^| or |$| \H\ asymptotically. 



i) Let us first consider the possibility ^ (x H. In that case ( p. 81) and (p. 91) imply that either 
\H\ ':$> or H (X H"^. In the first case we get i ~ 5H and then (|2?7|) imphes that 



-8^2 = 2K?oe^*e 



(4.2) 



i.e. negative energy. We thus exclude that possibility. In the second case from 
(|2.8fj2.9|) we obtain the solutions of the pre-big bang scenario, i.e. Eq. (|3.3|) . For 



consistency, we must verify that these solutions satisfy the constraint (|2.7D . This is 
the case only when the energy density is subdominant compared to the other terms 
in this equation. By using (^^), one can see that e^*e goes like -uzfi while the other 



terms behave like 



\t-to\ 

tQ. Thus the solutions 



jj—j-^, so the energy is subdominant as t 
can be approximated to those of the pre-big bang scenario near the pole, Iq. At the 
same time we see that this possibility is excluded for |t| — > oo. 



ii) For |6| < Eqs. (I]!) and (U) imply 27H'^ + 5H 
p = in Eq. (|2.11D gives us H + 9H^ 
excluded. 



0. On the other hand, setting 
0, a clear contradiction. Thus, |<I>| <C \H\ is 



iii) We are thus left with the sole possibility |<^| ^> \H\. In that case 



implies 



2$' 



2 2^ 

KiQe e, 



(4.3) 



and (U) together with (|2j|) imply + = 0, which gives $ 
With p = 0, Eq. ( ^.11]) then implies H - 2H^ = 0, which gives H - 



log(lt-tol)- 
where h 



{t-toV 



is some constant. This is consistent with Eq. (O) only for |t| 
as Itl — > oo. 



oo because e^*e 



1 



by 



To summarize, we find that the only consistent asymptotic solution for |t| — > oo is given 



$~-log(|t|) 



H ~ 



t 



2 • 



(4.4) 



and (U) for t ^ to. Note that i^E ^ -f for (g^). 

Now we are in the position to predict the possible string frame bounce scenarios. Fol- 
lowing the same logic as in the last part of the previous section concerning the NEC in the 
Einstein frame, one can show that the only possible bounce scenario is the transition either 
from pre-big bang-like solution'^ to pre-big bang solution or from post-big bang solution to 



Note that what we refer to as pre/post-big bang-like solution is given in Eq. (0) with negative/positive 
time 
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post-big bang-like solution. This is because the other transitions correspond to a bounce in 
Einstein frame, which is excluded in our model which satisfies the NEC. To explain this in 
more detail, we consider the case where we start with contracting pre-big bang-like phase. 
For large negative times our bounce solution is in agreement with the pre-big bang-like 
solution with accelerated contraction of the universe and growing dilaton. Then the uni- 
verse goes through a bounce and for t ^ to it approaches a pre-big bang solution with 
accelerated expansion and growing dilaton. The Hubble parameter in the Einstein frame 
remains negative and keeps decreasing. As will be seen in the next section, the numerical 
solutions are in good agreement with this picture. 



5. Numerical results 



In this section, we numerically solve Eqs. (^- 2.10| ) to obtain global solutions. In what 



follows we set 2\/2KlQTg = 1 (this can always be achieved by adding a suitable constant to 
the dilaton), so that 

2^?oe = e-*e-^'+^' (5.1) 
2^l^p = -e-*e-^' l^e^' - d^j . (5.2) 

With this setup, we performed the numerical analysis and found a family of bounce 
solutions. For example. Figs. § show a bounce solution with the initial conditions, ^>(0) = 
-5, 6(0) = 0.05, T(0) = 1000 and i?(0) = -0.002. The graph in Fig. |a shows the 
evolution of the Hubble parameter (solid line). The bounce takes place near t = 8. The 
bouncing solution can be seen as a transition from the contracting pre-big bang-like phase 
(short dashed line) to the expanding one (long dashed line). Both asymptotic solutions are 
obtained by setting p = in the equations of motions since the pressure is negligible in the 
far future and past (see Fig. ^b). The 'double bump' feature of the equation of state can 
be understood by noting that as the non-BPS brane builds up \T\ decreases and thus w 
increases from zero as explained in section 2. Then as T reaches the top of the potential 
|r| becomes small and consequently w decreases again. Indeed for T = the equation of 
state is that of a cosmological constant. 

We found that a broad range of the initial conditions are allowed for the bounce, so 
there is no fine-tuning problem. For instance, T(0) = ±10^ and T(0) = ±100 (keeping 
the same initial conditions as above for the other variables) gives bounce solutions with 
essentially the same behavior. Note that the large initial values for T do not represent a 
fine tuning. Rather it reflects the condition that the non-BPS brane is absent at very early 
times. We see that the asymptotic behavior of this family of solutions is similar to the 
expanding pre-big bang case, in which the Hubble parameter blows up. This agrees with 
the results of the previous section. 

If one changes the sign of $(0), a very different kind of bounce is obtained. For instance. 
Figs. ^ show a bounce solution with the initial conditions, $(0) = —5, 6(0) = —0.15, 
r(0) = 5 and i?(0) = -0.005. The graph in Fig. | shows the evolution of the Hubble 
parameter H . The bounce takes place near t = 2, and the universe smoothly evolves to the 
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^ t 



5 10 15 



(a) Hubble parameter H (solid line) (b) equation of state w 




(c) dilaton $ (d) tachyon T 



Figure 3: Bouncing numerical solution with the initial conditions $(0) = —5, $(0) = 0.05, 
T(0) = 1000 and H{0) = -0.002. 



standard cosmological regime, where the Hubble parameter and the dilaton decreases with 
time (Fig. I). But going back in time further, we found a singularity where the Hubble 
parameter blows up. This solution can be seen as a transition from the contracting post- 
big bang phase (short dashed line) to the expanding post-big bang-like phase (long dashed 
line). Both asymptotic solutions are obtained by setting p = in the equations of motion 
since the pressure is negligible in the far future and past (see Fig. ^). 

In addition, we found oscillatory solutions, in which a double bounce takes place (see 
Figs. P). The solution in Figs. |6| can be seen as a time reflected one of the solution of 
Figs. ^. Which solution is obtained depends on the sign of ^(0). In both cases the asymp- 
totic behaviors towards the curvature singularity are analogous to the non-oscillatory cases 
mentioned above. The evolutions of the equation of state are shown in Figs. I and |. Note 
that the negative equation of state implies small |r| (see Fig. ^), and this means that the 
oscillatory solutions can arise if the speed of the tachyon is small around the top of the 
tachyon potential (i.e. near T = 0). Alternatively this kind of solution can be obtained 
when we arrange |6| to be small enough around the top of the tachyon potential (see Fig. 0). 
The variation of <I>(0) with the other initial conditions fixed (in this case at $(0) = —5, 
H(0) = 0, and r(0) = 0) shows that decreasing <I>(0) gives rise to a transition from single 
bounce to cyclic bounce. 
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(c) dilaton $ (d) tachyon T 



Figure 4: Bouncing, numerical solution with non-singular future using initial conditions <&(0) — 
-5, $(0) = -0.15, r(0) = 5 and H{Q) = -0.005. 



H w 




(a) Hubble parameter H (b) equation of state w 



Figure 5: Oscillatory, numerical solution with initial conditions $(0) = —5, H{0) = 0, $(0) = 0.05 
and T(0) = 0. 

In conclusion, we obtained bounce solutions, where the universe smoothly evolves from 
the contracting phase to the expansion. The bounce scenario that we here found can be 
classified into the following two cases. Note that these are exactly in agreement with our 
predictions on the bounce scenarios in the previous section: 
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(a) Hubble parameter H 
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Figure 6: Oscillatory, numerical solution with initial conditions $(0) = —5, H{Q) = 0, $(0) 
-0.05 and T(0) 0. 
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Figure 7: Hubble parameter H for $(0) = 0.06 (solid line) and $(0) = 0.072 (dashed line). 



1. Transition from accelerating contraction (the contracting pre-big bang- like phase) to 
accelerating expansion (the pre-big bang inflation): In this case the dilaton grows up, 
and if the speed of the tachyon (or dilaton) is small enough near the maximum of 
the tachyon potential, a double bounce can take place (Figs. Isj) before the universe 
evolves to the pre-big bang phase. 

2. Transition from decelerating contraction (the contracting post-big bang phase) to 
decelerating expansion (post-big bang like phase): In this case the dilaton decays, and 
a double bounce can also take place (Figs. |^ under the same condition mentioned 
above. 

In all cases the tachyon rolls over the top of the potential in the course of its evolution, 
and the bounce seems to happen when the tachyon reaches around the top of the potential. 
The pressure is important only around the bounce, and negligible (dust) asymptotically. 

Our string frame bounce solutions correspond to monotonously contracting (i. e. = 
H — ^ < ) or expanding geometries {He > 0) in the Einstein frame (see Figs. |8[]9|), 
meaning that there is no bounce in the Einstein frame for our solutions. This is the fact 
that our model does not violate the NEC in the Einstein frame. 
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The numerical results presented in this section verify the results on the asymptotics 
in the previous section. Interestingly, as far as the dilaton $ is concerned the asymptotic 
solution obtained here agrees qualitatively with the global numerical solution. In other 
words the dynamics of the dilaton is not much affected by the presence of the non-BPS 
brane and is qualitatively similar to that of the pre-big bang scenario. Concerning H{t) 
things are different: See Fig. |3|for example. For large negative times H{t) is well described 
by the asymptotic solution described in Eq. (^^). Then near the bounce which takes place 
at t ~ the non-BPS brane affects H{t) significantly. Then for i — > to > (near the 
pole), where the pre big-bang singularity occurs, the Hubble constant is well described by 
the solution in the pre-big bang scenario. This can be understood from the fact that the 
brane has already decayed for t — > to- 



6. Nonsingular solutions - Example 

As we have seen in the previous section, there is a singularity either in the future or in the 
past, depending on the sign of and we expect that this singularity may be resolved in 
the same way as in pre-big bang scenarios, e. g. relying on a' corrections or quantum loop 



corrections or alternatively using a dilaton potential (see ||T2| for some explicit examples). 
As a matter of fact, resolving this kind of asymptotic singularities is less difficult than the 
big-bang singularity. In this section we will give an example of resolution of this asymptotic 
singularity. 

Near the singularity the curvature and the dilaton blow up, and this suggests that a 
potential term of the form -Re* in the Lagrangian may smooth out the singularity. Here 
the coupling between the Ricci scalar and the dilaton is introduced because the singularity 
appears both in the curvature and dilaton. Such a term is quite likely to appear as a' 
correction in the open string sector, since it has the form of a tree level correction in the 
open string coupling constant and R is the natural invariant built from background metric 
derivatives^ . 

■*In fact one can always arrange the bounce to take place at t = by shifting the time variable. 
'a similar potential has been motivated in the context of string gas cosmology in | p2| as a Casimir-type 
potential. 
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Thus as an example, in which such an additional term may resolve the singularities, 
we study the dynamics of the system where the action (|2.1|) is supplemented by a potential 



d^^x^/^e-^'^RVi<^). (6.1) 



The equations of motion then take the form 



72H^ (1 + y($)) + - 36H^ (^1 + - -^^^-^J - 2K?oe^*e = (6.2) 

2$ - 2^2 + 18H^ - (dH + A5H^^ ^1 + V{^) - ^^-^^^ " ^^106^*^/2 = (6.3) 

($ + 8H^) {V'{^) - 2V{<^) - 2) + 6^ (y"($) - 4y (^>) + 4y($) + 2) + 

(1 + V{^)) (36^2 + 8H) + K?oe^*p = (6.4) 
€ + 9H{e + p) -^p = 0. (6.5) 

We require that the additional term should not spoil the bounce, namely this term is 
important only when the curvature becomes very big. For concreteness we choose V = 
— e*^^/40. Using the same setup as in the previous section, we perform the numerical 
analysis. 

First, let us consider the case in which the dilaton grows; in this case we faced a future 
singularity. We found that the addition of a potential (|0]) can resolve the future curvature 



singularity. This is shown in Figs. 10, where the dashed (solid) curve corresponds to the 
case without (with) the additional term. Here the initial conditions are chosen such that the 
bounce takes place at t = (in other words we impose the initial conditions at the bounce 
and extrapolate in both directions in time). In the case without the additional term, there 
is a singularity, while in the other case the universe evolves to the standard cosmological 
regime where the Hubble parameter decreases. As can be seen from the plot, the dynamics 
are almost the same in both cases before the Hubble parameter gets significantly big, so 
that the bounce is not spoiled. Once the Hubble parameter is large enough, the additional 
term smooths out the singularity. 

Now we turn to the case in which the dilaton decreases. In the previous section we 
have seen that in this case there is a past singularity. With the help of the additional term 
mentioned above, we found that this singularity can be resolved as well. This is shown in 



Figs. [11|. The mechanism of resolving the singularity is analogous to the case of growing 
dilaton that we have seen above. What is interesting is that this solution corresponds to 
the time reflected version of the case of growing dilaton since H i— > —H and $ i-^ — $ under 
the time reflection, t —t. 

In both cases, the dilaton dynamics still has a singularity either in the future or in the 
past. This is in contrast to the case where the tachyon sector is absent (i.e. e = and 



p = 0). (see Figs. |12|, where the pre-big bang singularity is resolved not only in H, but 



also in <I>). Since the singular behavior of our solutions is the same as in the pre-big bang 



scenario, we expect that in principle the dilaton singularity can be resolved as in Figs. |12 
but this may require fine-tuning. 
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(a) Hubble parameter H (b) dilaton $ 



Figure 10: Numerical solution with growing dilaton (solid lines: with additional term, dashed 
lines: without additional term). 




(a) Hubble parameter H (b) dilaton <& 



Figure 11: Numerical solution with decreasing dilaton (solid Hues: with additional term, dashed 
lines: without additional term). 



To sum up, we have shown an explicit example in which an additional term that might 
arise from higher order corrections can resolve the curvature singularity without affecting 
the bounce dynamics, though the singularity in dilaton has not been resolved. 



7. Conclusion 



We suggested bounce scenarios, in which a non-BPS space-filling D9-brane in type IIA 
superstring theory drives a bounce of the scale factor in the string frame. We employed 
the lowest order effective action for the metric and dilaton in the string frame as well as an 
effective action for the open tachyonic mode of the non-BPS D-brane. The positivity of the 
pressure of the tachyon field is responsible for the bounce, which is why the DBI action, 
for instance, can not drive the bounce. The curvature as well as the time derivative of the 
dilaton remain small during the bounce. In other words, the gravitational sector is entirely 
classical. 
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Asymptotically our bounce solutions look like pre-big bang or post-big bang solutions, 
with singular behavior of the curvature and the dilaton. The asymptotic string frame cur- 
vature singularity can be resolved by adding a phenomenological potential, oc i?e~*, which 
may or may not result from a' corrections in the open string sector. It would be desirable to 
determine the sign and the precise numerical value of the proportionality coefficient. With 
our choice of the sign the gravitational coupling changes sign in the string frame. This re- 
sults in a bounce in the Einstein frame at some time after the bounce has taken place in the 
string frame without violating the null energy condition. An interesting observation is that 
while our phenomenological potential stabilizes the dilaton within the perturbative regime 
it fails to do so once the tachyonic sector is included. An obvious question is then whether 
a modified potential exists which stabilizes the dilaton in our model, and if so, whether it 
can be derived from string theory. We should also mention that throughout this paper we 
assumed the isotropy in 9-dimensional space (modulo compactification, the details of which 
did not play a role here) during the era of the bounce. For phenomenological reasons it 
may be preferable to consider scenarios with a different dynamics for the scale factor of the 
internal dimensions. In particular, orbifold compactifications are interesting since they are 
accompanied by a reduction of supersymmetry. A preliminary analysis shows that for a 
/Z2 orbifold the asymptotic solutions in the far past are not modified qualitatively. The 
numerical evolution of the global solution requires more work, however. We hope to report 
on this issue in a future publication. 
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